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^ ■ Abstract 

trr* I The field equations of the scalar-tensor theories of gravitation are presented in dif- 

■^ ' ferent representations, related to each other by confornial transformations of the metric. 

OS ■ One of the representations resembles the Jordan-Brans-Dicke theory, and is the start- 

O , ing point for the generation of exact electrostatic and magnetostatic exterior solutions. 



^- 



The corresponding solutions for each specific theory can be obtained by transforming 



^_ ■ back to the original canonical representation, and the conversions are given for the the- 

ories of Jordan-Brans-Dicke, Barker, Schwinger, and conformally invariant coupling. 
The electrostatic solutions represent the exterior metrics and fields of configurations 
/\^ • where the gravitational and electric equipotential surfaces have the same symmetry. A 

^ . particular family of electrostatic solutions is developed, which includes as special case 

the spherically symmetric solutions of the scalar-tensor theories. As expected, they 
reduce to the well-known Reissner-Nordstrom metric when the scalar field is set equal 
to a constant. The analysis of the Jordan-Brans-Dicke metric yields an upper bound 
for the mass-radius ratio of static stars, for a class of interior structures. 
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I Introduction and Summary 

The scalar-tensor theories of gravitation^^'^ differ from Einstein's general relativity since they 
admit the existence of a long range scalar field coupled to matter in such a way that the 
theories satisfy the weak principle of equivalence, and therefore are viable alternatives for 
the explanation of gravitational phenomena. However, the strong principle of equivalence is 
not valid, since the scalar field connects the local physics with the mass-energy of the rest 
of the universe. In fact, the Jordan-Brans-Dicke[JBD] theory, which has the simplest scalar 
coupling, was formulated as a relativistic model allowing a better representation of Mach's 
principle and Dirac's ideas concerning the connection between the Newtonian gravitational 
constant and the age of the universe^. The roots of the scalar-tensor theories go back to the 
original Kaluza-Klein five dimensional unification of the gravitational and electromagnetic 
forces, and the modifications by Einstein-Mayer^, Jordan^, and Thiry^. The five dimensional 
metric of these theories has five additional degrees of freedom, where four of them serve to 
represent the electromagnetic vector potential while the other behaves as a scalar field in a 
four-dimensional geometry. Thus we have that a scalar field appears as a natural outcome 
of the intriguing Kaluza-Klein unification and, furthermore, is the simplest long range field 
by means of which the matter distribution of the universe can affect local physics. The 
scalar-tensor formalism is able to incorporate this hypothetical but interesting interaction 
consistently and in agreement with experiments, and hence opens a broader theoretical 
framework to approach the study of gravitation. 

Aside from the physical and philosophical motivations for considering the existence of 



scalar fields, they could also provide heuristic representations of matter fields which, due to 
their mathematical simplicity, help us gain physical insights. For instance, it is well known 
that an isentropic fluid with pressure equal to energy density can be represented by a scalar 
field. Moreover, they could also be interpreted as large perturbations, serving in the study 
of the stability of general relativity backgrounds. 

There are several alternative representations of the scalar-tensor theories which are re- 
lated to the original canonical representation, Eqs. (2.1) and (2.2), by a conformal trans- 
formation of the metric, where the conformal factor is a function of the scalar field. For 
example, in the Einstein-scalar representation the new metric, g^^, and the original metric, 
Qf^^, are related by 

The field equations for 'g , given by Eqs. (2.6) and (2.7), are as in Einstein's theory but 
with an added energy- momentum tensor for the scalar field. A more general transformation 

is 

~ -t 

9fj.u ~ g^u -I 



where the new scalar field is an arbitrary function of 0, and we arrived at the field equations 
(2.12) and 2.13) which are remarkably similar to the originals. Alternatively, we can reexpress 
the Einstein-scalar representation using now instead of as the scalar field to obtain Eqs. 
(2.17) and (2.18). A particular choice of the function 0(0) makes the theory look similar 
to the JBD theory in the sense that the second order equations for the metric, Eqs. (2.12) 
and (2.17), take the same form as the corresponding one in the JBD theory, with a constant 



parameter, ujq, playing the role of the scalar field coupling constant. When uJq is chosen to 
be zero the theory can be written as a five dimensional Einstein's field equations in which 
the geometry is independent of the fifth dimension (Eqs. (2.22) - (2.28)). 

The representation 'g , 0, is suitable for the mathematical analysis, and therefore is the 
starting point in Section III for the generation of exact scalar-tensor electrostatic solutions, 
starting from given static vacuum metrics in general relativity. This mapping gives the 
exterior fields of configurations of arbitrary shape, for which the surfaces of constant elec- 
trostatic and gravitational potential coincide. They can be converted into new solutions of 
the original canonical representation, 0, g^^,, for each particular theory, and the necessary 
transformations are then given for the specific theories of JBD, Barker, Schwinger, and for 
conformally invariant coupling. Solutions representing the exterior gravitational, electric, 
and scalar fields of oblates or prolates configurations are presented explicitly in Section IV. 
These solutions reduce to the Reissner-Nordstrom metric when the scalar field is set equal 
to a constant and the configuration is taken to be spherically symmetric. 

The spherically symmetric solutions are study in Section V, where it is found that, in 
contrast to general relativity, in the scalar-tensor theories the analysis of the exterior solution 
alone is not sufficient to establish and upper bound for the mass-radius ratio of the static 
configurations, and it is necessary to consider the interior configurations as well. This is 
done in detail for the JBD theory, where an upper bound is determined for equations of 
state where the volume integral of the trace of the energy-momentum tensor is positive. 

In the last section, exact scalar-tensor magnetostatic solutions are given in terms of 



vacuum metrics in general relativity. Then an explicit family of solutions is worked out, that 
represents prolates or oblates configurations immersed in magnetic and scalar fields. These 
solutions are scalar-tensor theories generalizations of Melvin's magnetized metric of general 
relativity. 

II The Field Equations 

The scalar-tensor theories field equations, without a cosmological term, can be written in 
the following canonical form: 






0,M0,!^ - ^9f,u(t)''^<P,, 



+ 0;^. - ^^.0:«, (2.1; 



{2uj + 3)0:^ = ^^ [v^rl = SttGT: - 0'>,a^, (2.2) 



-9 
where, as usual, commas and semicolons mean partial and covariant derivatives, respectively, 

Gfj^u is the Einstein tensor, and T^jy is the energy momentum tensor for all fields excluding 

the scalar, 0, and the metric, g^i,, fields. We are essentially adopting the notation convention 

of Misner-Thorne and Wheeler^. It follows from Eqs. (2.1), (2.2) and the Bianchi identities 

that we still have, as it is in general relativity, the energy-momentum conservation: 

T^^^ = 0. (2.3) 

Thus, in this representation, the scalar field does not enter the equation of motion, and 
consequently free falling test particles move in geodesies of the metric g^^i,. With the simplest 
choice of coupling uj = const, we get the well known JBD field equations. Other examples 



of interesting theories can be defined by 



ijj 



h 



k2<P + h 



(2.4) 



where ki, k2 and k^ are constants. For instance, if ki = k2 = —k^ = 1, we have the theory due 
to Barker^'^, while the case /cs = was motivated by Schwinger^^. Another attractive model 
is the conformally invariant curvature coupling^, in which ki = 3,k2 = —1, and k^ = 1. 

An alternative representation of these theories is provided by the conformal transformation^^ 



Jfiu 'ryt^i'i 



(2.5) 



which gives rise to the Einstein's scalar form of the field equations: 



Gfiv — SttT^j^ 



2a; + 3 



3 

. + - 



1. 
T 



^,^.^,. - - V*^'"*^-" 
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-W 



,a 






where 



and 



$ 



T 

M^ 1 ' 



(2.6) 
(2.7) 



(2.8) 



(2.9) 



The bar over G^^, means that g^^, is replacing g^^, and now indices are raised with g^^^. 
A more general transformation is generated by letting 



9^^^y 



<p9,iu ^ V 





(2.10) 
(2.11) 



where is an arbitrary function of 0. The new field equations can be written in the form 



UJ 







2cj + 3 



(p,^^(|),u - 2fi'M^'^'°<^,< 



GT, 



57r- 



Q J,,o 






,0! 



a 



where 



T T 



a 



din(f) 



2Co + 3 



2lu + 3 



a^ 



Yet, if instead of g^v, the Einstein-scalar metric (^^^ is used, we will obtain 



G^, = 8nGT^, + (co 


3\ r 

+ 2). 


2cu + 3 
J-Ti 


y-g^'^] 


= 87r 



1. 
2' 



^,^.^,u - T^V*^'"*^. 



GT^ ~ ^~ du 
o- ' d^ 



where 



$ 



An interesting choice for the new scalar field is 



$ 



V2cJo 



/ V2uj + 3rf$ + Const. 

'0 + 3 J 



cjo = Const; 



(2.12) 
(2.13) 



(2.14) 

(2.15) 
(2.16) 



(2.17) 
(2.18) 



(2.19) 



(2.20) 
(2.21) 



since it then can be shown from Eqs. (2.15) and (2.16) that uj = cjq, and therefore the field 
equations take the JBD form, except for the factor - multiplying T^ and T^ in Eqs. (2.13) 



7 



and (2.18). Nevertheless, here ujq is an arbitrary parameter, which, of course, for the JBD 
theory can be set equal to u. 

The peculiar value uq = yields a representation that can be cast in the following 
five-dimensional language: 

V=^^, ; /i,z/ = {l-4}, (2.22) 

1^5= (2.23) 

L= 0' (2.24) 

GT^,,= S^ ; G= const., (2.25) 

GT^5= (2.26) 



5 5 1 ~ 

GT55= i^4> 



1-i 

a. 



Gf^, (2.27) 



and hence the scalar tensor field equations become 



GAB=87rGTAB; A,B = {l-5}, (2.28) 

5 5 

where Gab is the five- dimensional Einstein's tensor for the metric Qab [X'^). 

We can readily integrate Eq. (2.20) for the u given by Eq. (2.4), and we obtain: 
(i) Barker's Theory 



2 

(ii) Schwinger's Theory 



-Tan-^^cj) - 1 + Const. ; £ = ^/2uq + 3, (2.29) 

f 
(f) = Tan^-[ln(t) - Const.] + 1. (2.30) 



2 k 

+ Const. ; a = -^ (2.31) 



(i^/axj) ki 



4 
-— [Const. -irKpr^. (2.32) 



(iii) Conformally Invariant Curvature Coupling 

£n(f) = —T-d-n 



1 + r 
4Const' 



Const., r= ^1-0 (2.33) 

(2.34) 



0VV3+ Const.'] ^ 

III Generation of Electrostatic Solutions 

Let's consider the field equations (2.17) and (2.18) in the representation where Co = ooq 
[JBD hke representation], and assume that we have the energy-momentum tensor of an 
electromagnetic field: 



InGT^, = 2 



a TT I?2 



F F" a F 



(3.1) 



4^ 
_dA^_dA^ 

Furthermore, let's assume that the space-time is static, and the vector potential only has the 
component A^. Hence, we are dealing with three potentials; the electrostatic potential A4, 
and the gravitational potentials which are represented by g^^ and 0. We know that when the 
configurations are spherically symmetric these potentials depend only on the radius variable 
that characterize the equipotential surfaces. Moreover, even if the configurations are not 
spherically symmetric, we could still expect that the gravitational and electric equipotentials 
are surfaces of the same shape, and, consequently, these potentials once more would be 
identified by a single function, say u. Thus, in searching for new solutions, it is quite natural 



to start with the following ansatz: 

Ai = A^iu), (3.3) 

= 0(«), (3.4) 

gu = 9uiu). (3.5) 

Then one finds (Appendix A) that the electrostatic solutions of the JBD-like field equa- 
tions are expressed in terms of solutions, g , of Einstein's static vacuum field equations, 
G = 0, in the following way: 

'''~[b-i-g'.,yr ^ ^ ^ 

g^, = ^^^^^^^^i-g'uY-%;hj = 1,2,3, (3.7) 

, _ ±Mb + (-^44)1 I , .^o^ 

$ = e £n(-^44) + Const., (3.9) 



where a, b, c, and d are constants. 

For asymptotically Minkowski space-time where 

2CM' I 

- ^L ^ 1 - —^ ; r = ^X,X^ ^ oo, (3.11) 

we have that, with the choice 

{b 

a ^ —- 

46 

-ii 
'[1-6]' 



-^, (3.12) 

rf = ±v^tt|, (3.13) 
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the asymptotic behavior of g , and A4 are 



- t/44 ~ 1 



2GM 



Uij ^~~' iJij 



■ 2GM' , 

1 + c-1 



2c 
6-1 



const. 



A, 



1- 



2G'eM' 



A;g 



(3.14) 

(3.15) 
(3.16) 
(3.17) 



where /c is a constant, and 



M = cM' 
kQ = 



b+l \ 
2cVbGM' 



6-1 

GM^ = {GcM'f + {kQf 



(3.18) 

(3.19) 
(3.20) 



Note that if 6 ^ 00, then Q -^ 0, and M = cM' . Thus for M and M' to be positive we must 
assume c > 0. From Eqs. (3.14) and (3.17) we see that M and Q can be identified with the 
total inertial mass and electric charge, respectively. 

The solutions (3.6) - (3.9) can be converted into the corresponding solutions of the original 
scalar-tensor field equations (2.1) and (2.2) by using the transformation (2.10): 



9fiv 



9 



flU 



(3.21) 



and then 0(0) could be determined by the inversion of Eq. (2.20), with given, by Eq. (3.9), 
explicitly as a function of (744. For instance, for the JBD theory we can have = 0, while for 
Barker's Schwinger's and conformally invariant curvature coupling theories, the necessary 
inversions are found in Eqs. (2.30), (2. 32), and (2.34), respectively. 
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IV Weyl Type Solutions 



An interesting family of solutions are generated with the transformation (3.6)-(3.9) when we 
take as seed metric, (7' , the following well known Weyl solution: 

r o/qT^ 

(4.1) 



(4.2) 



^44 


= 


1- 


2/31 
r 


S 




1 - 


2^1 


52-<5-l 






r 








1 _ M + ^stn^^ 



<52-l 



9'^^ = r^ 



1- 



2p- 



l-<5 



sin 9 



(4.3) 



(4.4) 



where j3 and 5 are constants. Clearly, the above solutions reduce to the Schwarzschild metric 
when 5 = 1. Otherwise, they represent the gravitational field exterior to oblates (5 > 1) or 
prelates (5 < 1) ellipsoidal configurations^^. 



The new solutions generated by Eq. (3.6)-(3.9) are: 

2 



6-1 



-^44 



1 



2/3 



5c 



r 



«■ 



r 



l-W. 



52_5c-l 



i/ee 



^, 



i, 



^[1 - f + frsm^^] 
r2[l 



M + 4eV^2fll52-l 

2[i 2/3]l-5c 

2/ 



A4 = ± 



6-l/6+[l-^ 



sin ^ 



2/3] 5c ^ 
r 



2V6 V&-[i-¥P 

12 



6+1 
2v^ 



(4.5) 

(4.6) 

(4.7) 

(4.8) 
(4.9) 



$ = (5e In 



where 



i, 



r 



6-1 



+ Const. 



(4.10) 



(4.11) 



yb-[i-fY-J 

The solutions (4.5)-(4.10) represent static gravitational, electric, and scalar fields with oblate 
or prolate equipotential surfaces (Appendix B). When 6 = 1 the solutions are spherically 
symmetric, and the metric takes this simple form: 



sin^O 



'h-[l 



2/3 1 



6-1 



1- 



2(3 



l-c 



r^ = R^ 



9aA ~ 9rr 



(4.12) 



(4.13) 



The above solution is a scalar-tensor generalization of the Reissner-Nordstrom metric, to 
which it reduces when c = 1 (Appendix C). 

If c 7^ 1, and r = 2f3 = 2GM', we see from Eq. (4.13) that g'^^ = 0. However, in contrast 
to general relativity where r = 2GM' represents an event horizon, here we have, from Eq. 
(4.12), that 

E-^0, 



lim r ^ 2(3 



(4.14) 



and, since the proper areas of constant radius of the new solutions are equal to A-kR , then we 
see that the Schwarzschild horizon is mapped into the origin of the JBD like representation. 
Thus, we have that the static metric (4.12), (4.13) and scalar field (4.10) describe the whole 
space < -R < oo, for 2GM' < r < oo, without a horizon. Furthermore, this exterior metric 
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does not imply a lower bound for i?, analogous to r > 2GM' appearing in general relativity. 
Nevertheless, we have to be cautious with the physical interpretation of R, since now we 
are working in the JBD like representation where the physical meaning of the metric and 
scalar field are not transparent, for we have that free falling test particles do not move in 
geodesies of the metric y . Instead, we should look at the behavior of the radius, R{r), in 
the original canonical representation, g^iy, 0, which, by virtue of Eq. (2.5), is related to R in 
the following way 

R = .fWo = M=^. (4.15) 



; — 9ee _ R 
V (p V0 



In the next section we will take on this analysis for the JBD theory. 



V Analysis Of R(r). 



For the JBD theory we can have cj = cjq? ' 



and hence from Eq. (4.10) we find that 



1 - 2GM' 



.6 = 1, 



(5.1) 



where the constant of integration has been chosen such that ^ 1 as r ^ oo. Using Eqs. 



(4.12), (4.15) and (5.1), we get 



R = r 


'b- 


Q 2GM'\c 




6-1 



2GM' 



(5.2) 



and depending on the sign of S 



1— c— e 



we will have as r ^ 2GM' that either i? -^ 0(5" > 0), 



or R ^ oo{S < 0) after reaching a minimum value for some r > 2GM'. On the other hand, 
the sign of S is determined by the value of c, which, in turn, depends on the interior structure 
of the star, as described below. 
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The field equations (2.17) and (2.18) with uj = ujq imply for a static system that 



-g[tn - ^44 j 



-g^)' 



-16nG 
8nGT 



7t=;4 _ Itt^q 

U 2^ a 



:« a{2uJo + 3) 



-9: 



(5.3) 
(5.4) 



and consequently 






d'S, = 8nG 



Tt - -t: 



nv^d^x 



-gd^X. 



(5.5) 
(5.6) 



(t(2cjo + 3) ■ 

Thus, if the two dimensional surface of integration, d'^Si, is chosen outside the boundary of 
the star, we can then substitute, in the left hand side of Eqs. (5.5) and (5.6), for ^44 and 
$ the expressions given in Eqs. (3.6) and (3.9). At this point we simplified the problem by 
assuming that the configuration is neutral {Q = 0), and consequently Eqs. (5.5) and (5.6) 
become 



c / V^i^n - g'^^Yd^Si = -16nG 

e f V^{in - g'^^yd^S^ = SnG f 
From which it follows that 






"^^3 



gd'X, 



I 



a(2a;o+3) 



a(2cuo + 3) ■ 



J a(2wo+3) 



(5.7) 
(5.8) 



-2/ 



't^ _ It" 



W^x -2/ 



T4 _ IT" 
^4 2^a 



^#X' 



(5.9) 



where we used Eqs. (2.5) and(2.9). For a perfect fluid we obtain 

r (3p-pU/^gd^ 
e _ J o-(2a;o+3) 



c J{3p + p),/^d^X- 



(5.10) 
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In particular, if the pressure, p, is proportional to the energy density, p = ep, we find 



(36 - 1) 



c {3e + l)^2uJo + 5\^2uj + 3 ^ 



(5.11) 



where we use Eq. (2.16) with uj = ujq, and 



: r- r p^f^d^X 

V2l^T3y ~ Jpy/^d^X' 



(5.12) 



Hence for the JBD theory we get, with Uo = 00, 



3e-l 



(2cj + 3)(3e + l) 



(5.13) 



from which, using Eq. (3.10), 



2uj + 3 



,2- 



(5.14) 



Y^ + ^ + iM)' 

If e is not a constant, or we do not have a perfect fiuid we can still use Eqs. (5.13) and (5.14) 
by replacing e in them by 



JTiy^d^X 3 



JTt^/^d^X 



+ 1 



(5.15) 



Using Eq. (5.13) and (5.14), we can express 5 as a function of uj and e: 



S = -{l-c-e) 



( 



1 



I ^ 



2cj + 3 



2cc; + 3 



3?-l 
3?+l 



1 + 



3e-l 



.\ 



(2cu + 3)(3e+l) 



(5.16) 



Thus, assuming that 2cj + 3>0, ifO<e< |, then S* > 0, and when e > |, one finds that 
S* < 0. The case e = | yields c = 1, which implies that = const, and we are back to the 
Schwarzschild metric of general relativity. According to Eq. (5.15), for a positive energy 
density, -T^ > 0, we have e > |(< |) ii JT^^/^d^x > 0(< 0). 

16 



Let us first consider the case e > |, S* < 0, where it can be shown that R has a minimum 
value at 

R = Rrmn = 2GM'^ — ^ ■ u = S > 0. (5.17) 

Consequently, the solution does not cover the totality of space, since it leaves out the region 

< R < Rmin- Then, in order to be able to match an interior static solution to this exterior 

solution, it must be required that the radius of the interior configuration is larger than Rmin, 

thus implying that 

R > Rrmn = ^^^(^ + ^)'^^ > 2GM, (5.18) 

c U'^ 

where we used M = cM'. The value of U can be calculated as a function of u for a given 
interior model using Eqs. (5.15) and (5.16). 

We now turn to the case e < |, S* > 0, which implies that 

R>0; r> 2GM', (5.19) 

and therefore no lower bound appears for R from the above considerations of the exterior 
space-time alone, and further interior analysis is necessary. 

Besides their mathematical interest and heuristic value, the solutions (3.6)-(3.9) provide 
a necessary tool for confronting with reality the predictions of a very large class of viable 
theories, since with the more refined observations of strong fields it is increasingly impor- 
tant to use exact results. In particular, the solutions of this paper serve to investigate the 
trajectories of photons and particles in the space exterior to static systems. For instance, if 
the configurations are spherically symmetric, then, from 

f7"f7„ = F'^fJ^ + TlVr + U^U4 = -1, (5.20) 
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F" = — ^; dr^ = -g^JX^'dX^ = <pdT^, (5.21) 

dr ^ 

we derived the following expressions for the orbits: 

, . , . r L dr , , 

t{r)-t{ro) = ± _ _ , , (5.23) 

where, without loss of generality, we have assumed that the motions occur in the plane 

9 = n/2, and 

L^l^l (5.24) 

U 4 

is the so-called angular momentum per mass or impact parameter. Since the field are inde- 
pendent of t and (/9, then the energy-momentum conservation, T'^"; z/ = 0, yields 

mUup 
P^ = mU^ = — — = const.] m = const., (5.25) 

02 

Tflll A 

Pi = mll^ + eAi = — J — h eA^ = const., (5.26) 



from which we obtain that 



L - , , 

=; L = const., (5-27) 



1 eAjL 

and thus, for a neutral particle L = L, and U4 = const. (p 2 . Therefore, with the metric (4.12), 
(4.13) and the expressions (5.22), (5.23), we are in a position to extract the values of ^p{r) 
and t(r), from computer tabulations, and hence to obtain precise predictions for relativistic 
celestial mechanics effects, like time delay, light deflection, perihelium shift, etc. 



VI Magnetostatic Solutions 

Let us now assume that the metric and fields are static and independent of a coordinate X"-. 
Then we can generate magnetostatic JBD hke solutions starting from an Einstein's vacuum 



static metric, (7' , in the following way (Appendix A): 



^44 



-^44)^ 



^d'YYiPri 



bm + 9aai-9u 



I \l-c 



2 _ (-^44)^ 



^71 



9aa = i^mi-944) ^9aa^ 



9ii = ; 9ii ; ^,J 7^ 4,a, 



A 



A -±,r- Waai-9'ur'-b^] 



9aai-9MV ' + &m 



$ = e ^n(— (744) + const., 



where 



i^r. 



^Ctrr), Of] 



(6.1) 

(6.2) 
(6.3) 
(6.4) 
(6.5) 

(6.6) 



In order for g^^, to reduce to g'^^, when c = 1 and Aa = 0, the constants a^, bm and dm should 
be related as follows: 



^m ^\(^my 



^m '^(^m- 



(6.7) 



(6.^ 



We can get explicit magnetostatic solutions by choosing as background the axially sym- 
metric Weyl metric (4.1)-(4.4). In this case, X" can be identified with the azymuthal angle 
ip, and we obtain 



^44 = -A' 



1 - 



2P' 



Sc 



(6.9) 



19 



A2[1_M' 



2/3l52-5c-l 



9ee 



9 
A 



w 



[1_M 


+ 72 


sin^ef'-^ 


AV2 


1- 


r ■ 


[1_M 

L r 




sm2^]^'-i 


r2[l- 


_ Ml 

r - 


l-(5c 




A2 




Er2[l 


_ M 

r 


l-<5c 

sin^O. 



A 

^ = 5e in 



r 



(6.10) 

(6.11) 

(6.12) 
(6.13) 
(6.14) 



where 



A = B'r^ 



r 



l-Sc 



sm 



B 



±1 
2,/a: 



(6.15) 
(6.16) 



When (5=1, c=l, and /3=0, the above solutions reduce to the well known magnetized 
Melvin solutions^^, which represents a cylindrically symmetric bundle of pure magnetic flux. 
If /3 7^ we have a Schwarzschild black hole embedded in the magnetic flux^^. If, furthermore, 
(5 7^ 1 and c 7^ 1, the conflguration is not spherically symmetric, and a scalar fleld is also 
present. 
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Appendix A 

Suppose that the JBD hke field equations are independent of the coordinate X", and that 
the metric and electromagnetic vector have the form 

V ^ ( 'o I, ) • (*■!) 

Then, with the following transformation of the metric 

9aa = a'aa^^^ (A.3) 

5« - % (A'^) 



we obtain, from the JBD like equations, 

G., - SttGT,, - (.;o + h[^,^, - ^^,,$''=$,d = G^ - ^[A'' + i?'1 = 0, (A.5) 



\)[^.^. - \^^,^''~^A = G[^ - ^^^^' - «"=! 



where 

1 



Cijik = 9ii9jk - 7;9ij9ik: (A.6) 



A'' = iingiyiini^y' + (£nC)''(^nV^)'' + {ini^Y (ini;)'' + -^^, (A.7) 

B^^ = {ing'^y{inny^ + {ing'^y^innY + {inny{inVLy^ + (2cuo + 3)l>•'<l■^ (A.8) 

C = 9aa^- (A.9) 

Furthermore, 

2 V-5'' fl'aa 
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-ipQ 


V^A, 


,fc" 


;^ 


V-9' 


[ 9'Li' J 






7-^'$'' 


,fc 



0, 



(A.ll) 
(A.12) 



Note that in the right hand side of the equations the indexes are raised with g''^K 



We are looking for functions tp and Vt such that A = 0, B = 0, and, thereby, G^ • = 0. 
The inspection of the expressions (A. 7) and (A. 8) suggest that we work with the following 



ansatz: 



Aa = Aaig'aa) , ^ = ^^(^la)> 



$ = ^gU ; fi = ^{g'aa)- 



(A.13) 



(A. 14) 



The hipothesis (A.13) and (A. 14) are equivalent to the assumptions (3.3)-(3.5), and hence 
these are solutions where the surfaces u = g^^^ = const, are equipotentials of Aa^g^a and 
0. Thus, for instance, for static systems, a = 4, the surfaces of constant gravitational and 
electric potentials coincide, which is to be expected in astrophysical systems in equilibrium. 
Aside from the assumptions (A.13) and (A. 14) the functions ip and ^2 are arbitrary and 
can be chosen to satisfy the following relations 



dim/j 
d^ng'^a 



-l±v/l-4C#?, 



where 



diriQ 
ding'^^ 



kl 



-1 ± xll - e 



2-1 



(A.15) 
(A.16) 



1 dAg 

i' dg'aa 



,9 / . f d^ 

kl = {2uJo + 3) I - 



\ding'^^ 



(A.17) 
(A.18) 
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which can be shown to guarantee that A^'' and B^'' vanish. Then (A. 5) imphes that Gj • = 0, 
and as a consequence of the Bianchi identities, G.'^^'^ = 0, we also have G^^ = 0, or G = 0. 
On the other hand, from Einstein's vacuum equations, G = 0, we get 

1 



R = 



Jaa) 



= 0, (A.19) 

k 



which together with Eqs. (A. 11), (A.12), (A. 16) and (A. 17) imply, if ^^ > 0, that 

ki = const., (A.20) 

k2 = const. (A.21) 

Then from Eqs. (A.15)-(A.18), we obtain 

9aa = g'L^ = ibfgT^^ ' ^' ^ i' ^ = "^'''^•' ^^-^^^ 

Aa = ±^^§tL^ + const, (A.23) 

[b + 9aa] 



n = g',,^-\c^±^l-kl (A.24) 

V2cJo + 3 
We can also verify that Eq. (A. 10) is satisfied by the above equations. 

Equations (A.22)-(A.25) and Eq. (A. 4) can be converted into Eqs. (3.6)-(3.9) by identi- 
fying 'a' with '4'. If, on the other hand, x" is a space-like coordinate, then Aa represents a 
magnetic field. 

If the fields are also independent of another coordinate x'', then we can consider instead 
the conformal transformation 

9aa = %^, (A.26) 
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^.. = ^, (A.27) 

^., = ||;^,jVa,&- (A.28) 



and then Eqs. (A. 5) are derived with the following changes in A''-'' and 5^^; now in A^^, g\ 



means 

9'L = ^; (A.29) 

and in B^^, g'^^ is substituted by g'^^^^. Consequently, following the same routine that led to 
Eqs. (A. 22) - (A.25), we will again obtain Eqs. (A.22) and (A. 23), but now 

n = g',r\ (A.30) 



and 

£nl> = ^ Ing'uu + const. (A.31) 

Choosing 6 = 4 gives rise to the magnetostatic solutions of Section VI. If, on the other hand, 
Xb is space-like, then we would get electrostatic (a = 4) or magnetic solutions (a ^ 4) that 
were not studied in this paper. 
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Appendix B 

We see from Eqs. (4.5), (4.9) and (4.10) that the surfaces with r = const, are equipotentials 
for the gravitational, electric, and scalar fields. It is easy to calculate the proper distance on 
these surfaces along the circumferences 9 = const., and we get, using Eq. (4.8) 



l-(5c 



C^[e) = / x^Zdip = ^ — ^ sine = 27iRsine. (B.l) 



On the other hand, the proper distance around the curves ip = const, are given by the 
integral 



p^sin^e ^ ' '' 



Ce = 2 I ^QeedO = 2R I ( 1 + ^^ ^ dO (B.2) 



/o * JQ \ r'^[l 

where we used Eq. (4.7). We can see that since r > 2/3, 

Ce > C^{90°) = 2nR ; < 6 < 1, (B.3) 

Ce< 0^(90°) , (5>1, (B.4) 

and therefore the equipotential surfaces have shapes like prolates (0 < 5 < 1) or oblates 
{6>1). 
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Appendix C. The Reissner-Nordstrom Solution. 

If c = 1, one finds from Eq. (4.12) that 



r = :R-^, (C.i: 



and then Eq. (4.13) becomes 



-5«^f-^l-^ + ^^ (C,2) 

On the other hand, recalhng Eqs. (3.18) and (3.19), and since c= l,/3 = GM', we find that 

_ _, 2GM Q' 

- 9U = 9rr =1 ^— + «^ =2 ' (C-3) 



which together with Eq. (4.12) is the standard form for the Reissner-Nordstrom metric. 
Furthermore, Eq. (4.9) with 6 = c = 1 gives the famihar electrostatic potential: 

A, = ^ (C.4) 
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